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Final-state-exclusive momentum distributions of fast, forward travelling residual nuclei, following 
two nucleon removal from fast secondary radioactive beams of projectile nuclei, can and have now 
been measured. Assuming that the most important reaction mechanism is the sudden direct removal 
of a pair of nucleons from a set of relatively simple, active shell-model orbital configurations, such 
distributions were predicted to depend strongly on the total angular momentum / carried by the 
two nucleons - the final state spin for spin + projectiles. The sensitivity of these now- accessible 
observables to specific details of the (correlated) two-nucleon wave functions is of importance. We 
clarify that it is the total orbital angular momentum L of the two nucleons that is the primary 
factor in determining the shapes and widths of the calculated momentum distributions. It follows 
that, with accurate measurements, this dependence upon the L make-up of the two-nucleon wave 
functions could be used to assess the accuracy of (shell- or many-body) model predictions of these 
two-nucleon configurations. By use of several tailored examples, with specific combinations of active 
two-nucleon orbitals, we demonstrate that more subtle structure aspects may be observed, allowing 
such reactions to probe and/or confirm the details of theoretical model wave functions. 

PACS numbers: 24.50.+g,23.20.Lv, 21.60.Cs 



I. INTRODUCTION 

The momentum distributions of the residual nuclei, fol- 
lowing the removal of a single nucleon from a fast ra- 
dioactive secondary beam, offer sensitive probes of both 
strongly-bound and weakly-bound single-particle struc- 
ture near the (asymmetric) Fermi surfaces of neutron-rich 
and neutron-deficient nuclei. Specifically, the shapes and 
widths of the exclusive residue momentum distributions 
were shown to be characteristic of the orbital angular 
momentum of the removed nucleon [J-Q- 

The simplest generalization to the case of direct two- 
nucleon removal is to describe the wave function of the 
two nucleons in the projectile by a product of nucleon 
wave functions in assumed single-particle orbitals. Do- 
ing so, the two nucleons are uncorrelated, other than both 
being bound to the same core d, 0] ■ The heavy residue 
longitudinal momentum distributions in this limit, being 
essentially the convolution of those of the single-nucleons, 
depend on the assumed quantum numbers of the two nu- 
cleons, but, in the absence of explicit antisymmctrization 
or total angular momentum coupling of the two nucleons, 
are not characteristic of specific residue final states Q ■ 

More recent theoretical developments now treat fully 
the shell-model correlations of the two removed nucle- 
ons in the projectile many-body wave function [H, H[ . In 
the fully-correlated models the product of nuclconic wave 
functions is replaced by the shell-model two-nucleon over- 
lap, incorporating (i) the two-nucleon parentage coeffi- 
cients with respect to each residue final state (the two nu- 
cleon amplitudes, or TNA), (ii) proper antisymmctriza- 
tion of the two removed nucleons, and (iii) proper angular 
momentum coupling. 

The resulting theoretical description, and the insights 
developed here, are equally valid for reactions that re- 



move two loosely- or strongly- bound nucleons. How- 
ever, as has been discussed elsewhere @, in the case of 
removal of two loosely-bound nucleons the direct removal 
cross sections, of interest here, will be overwhelmed (ex- 
perimentally) by indirect reaction (one-nucleon removal 
plus evaporation) events. See for example reference (lOj 
for a quantitative consideration of the direct and indi- 
rect two-neutron removal reaction contributions along 
the neutron-rich carbon isotopic chain. For these rea- 
sons we will restrict our attention to examples for which 
the removed nucleons are strongly-bound, where the in- 
direct removal paths are effectively closed, and for which 
the direct cross sections are accessible experimentally. 

Demonstrative test cases, e.g. in Ref. pd| . assumed 
the two nucleons originated from a single orbital - a pure 
configuration. In this limit the TNA enter only as a mul- 
tiplicative (spectroscopic-like) factor and thus the new 
and interesting characteristics of the residue momentum 
distribution are a result of the correlations due to anti- 
symmctrization and angular momentum coupling. 

These developments demonstrate the potential of two- 
nucleon removal for exotic nucleus spectroscopy, show- 
ing the final-state-exclusive residue nucleus momentum 
distributions to have shapes and widths that are char- 
acteristic of the total angular momentum, /, carried by 
the removed pair of nucleons - and permitting final state 
spin assignments to be made [l2[. For the spin Ji = + 
projectile nucleus examples used in Ref. [ll| , there was a 
high sensitivity of the residue momentum distributions to 
the final state total angular momentum Jj = I. More- 
over, the shapes of these calculated distributions were 
robust with respect to variations of other key structure 
and reaction parameters, such as the nucleon separation 
energy. Further, the consideration of pure configuration 
examples (e.g. two protons, assumed removed from a 
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single active 7rld 5 / 2 , 7rld3/2 or 7r2si/2 orbital) showed 
considerable insensitivity of the two-nucleon removal dis- 
tributions to these individual nucleon quantum numbers; 
in stark contrast to results from single-nucleon removal 
reactions where the orbital angular momentum is critical. 

Thus, although the two-nucleon removal process is 
powerful for final-state spin spectroscopy in very exotic 
systems, its sensitivity to and ability to probe finer de- 
tails of the shell-model wave functions and the two nu- 
cleon configurations and correlations therein remains less 
clear. Already in Fig. 5 of Ref. Q, for the case of 
two-proton removal from 28 Mg, using the fully-correlated 
shell-model wave functions for the transitions to the first 
two low-lying 26 Nc(2 + ) final states one observed momen- 
tum distributions with different widths; demonstrating a 
sensitivity beyond the final state spin. Our objective here 
is to elucidate this sensitivity of the calculated residue 
momentum distributions to the particular two-nucleon 
configurations present and to understand the sensitivity 
to the combination of orbitals involved for a given value 
of the pair's total angular momentum /. 

In contrast to two-nucleon transfer reactions, such as 
the (p, 3 H) reaction, wherein the (p| 3 H) light-ion struc- 
ture vertex preferentially selects the pick-up of a spin- 
singlet (S ~ 0) neutron pair, the two-nucleon removal 
mechanism is not explicitly selective in the nucleon spins 
Q . Both spin-singlet and spin-triplet components of the 
two-nucleon overlap will be probed and, under the as- 
sumption that the residue and nuclcon-targct interac- 
tions (and S'-matrices) are spin-independent, the S = 
and S = 1 terms contribute incoherently to the reaction 
yield. We will show specifically, in the same way that 
single-nucleon removal is sensitive principally to the or- 
bital angular momentum £ rather than the total angular 
momentum j of the nucleon, that the two-nucleon re- 
moval reaction momentum distributions are sensitive to 
the components in the two-nucleon overlap with a given 
value of total orbital angular momentum L, L = £i 
The presence and relative strengths of these L compo- 
nents are determined via the shell-model overlaps and 
their TNA. Since, for a spin-zero projectile, the residue 
has total spin J/, with J f = L + S, the total spin con- 
tent of the overlap is determined by the nuclear structure. 
However, all spin components present are sampled by the 
reaction mechanism. 

Recognizing the sensitivity to L allows greater probing 
of the shell-model wave function, particularly in states 
where the mixing of several available nucleon configura- 
tions may be weak. Specific examples can exhibit par- 
ticularly strong sensitivity to the orbital combinations 
of the pair. In cases of strong mixing and sharing of 
strength among several orbitals a particular value of L 
may nevertheless be favored. In such cases the residue 
momentum distribution will be characteristic of the de- 
tails of the states populated and not simply of the total 
angular momentum value / of the nucleon pair. Of par- 
ticular interest will be different regions of A and Z that 
affect the originating, active orbitals of the two nucleons. 



In Section [IT] wc outline the approximations assumed 
and then present the required formalism for residue mo- 
mentum distributions within the LS'-coupling scheme. 
We will retain isospin labels for clarification of the un- 
derlying symmetries. In Section IIIII the importance of 
the total orbital angular momentum is elucidated by a 
detailed consideration of the spatial and angular corre- 
lations of the two nucleons that arc inherent in the two- 
nucleon overlap function. This analysis also provides in- 
sight into the possible sensitivity to the mixing of orbitals 
across major shells. In Section HVl we then consider par- 
ticular examples with projectiles of different A and Z 
where interesting effects are predicted. Examples will 
look at specific final states that could be populated in 
two-nucleon removal from the p-shcll, 12 C(— np), the sd- 
shcll, 26 Si(— 2n), and also the sdpf-cross shell situation 
in 54 Ti(— 2p). Wc summarize the article and draw con- 
clusions in the final Section. 



II. FORMALISM 

We discuss the sudden, direct removal of two nucleons 
from a fast projectile beam incident on a light nuclear tar- 
get at energies of order 80 MeV per nucleon and greater. 
In this intermediate energy range there have been exten- 
sive (positive) assessments of the validity and accuracy 
of the sudden/adiabatic and eikonal reaction dynamical 
approximations (see e.g. Section 3.5 of Ref. [l3[ and ref- 
erences therein) and of the theoretical ingredients used in 
their implementation; such as the importance of Pauli- 
blocking on the effective nuclcon-nuclcon interaction used 
141 ] . These wide-ranging assessments, carried out within 
the one-nucleon stripping and breakup contexts, remain 
valid in the present analysis; e.g. the role of strong ab- 
sorption between the projectile and target in reducing 
the effective reaction time and the energy for validity of 
the sudden/fast adiabatic approximation (TBI ]. 

We first briefly discuss the salient features of the ap- 
proach, previously developed in detail in Rcfs. d, |^, □J. 
Our emphasis here will be on the use of the LS'-coupling 
representation to describe the two-nucleon structure 
overlaps and to derive the expressions for the residue mo- 
mentum distributions in this basis. We will consider only 
those (stripping) reaction mechanism events in which the 
two removed nucleons interact inelastically with the tar- 
get nucleus. The role on the cross sections and mo- 
mentum distributions of the other major class of events 
(diffraction-stripping), where one of the nucleons is re- 
moved by an elastic interaction with the target, were 
discussed fully in Refs. @ and [llj, respectively. As 
was shown there, events from this second mechanism give 
residue momenta that are essentially identical to those of 
the stripping mechanism. These conclusions remain un- 
changed and will not be repeated here. 

We assume the projectile nucleus to be an antisym- 
metrized A + 2-body (shell-model) system, denoted by 
= ^jtMiTiTi (A, 1, 2), carrying total angular momen- 
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turn Ji and isospin Tj. In a high-speed collision with the 
light target, two nucleons may be removed to produce 
an A-body reaction residue in a final state /, often re- 
ferred to as the core state for simplicity. This final core 
state is & F > = §j f M f T f Tt {A). Each residue final state 
is denoted by /, while F = (/, Mj) is used to refer to a 
state with a specific angular momentum projection Mf. 
One is reminded that the final state of the two nucleons 
and of the target nucleus are unobserved and that the 
observables discussed are inclusive with respect to these 
degrees of freedom. 

A. Two-nucleon overlap 

The direct reaction will probe the two-nucleon overlap 

M F) ^1^(1,2) 

= (^(A)\%(A,l,2)) 
= E C?{IpJ f M f \JiMi) 

IfiTa 

(TrTfTflTin) [^(1)^(2)17; , (1) 



with D a = 1/ y/2(l + 5/3 1 j3 2 ) and where the angular mo- 
mentum and isospin couplings used are summarized in 
Fig. [TJ The nucleon-wave functions tpp{ r i) are 

^(r i )=Mn)Y iX (f i ) ■ (3) 

It is convenient to combine the statistical factors and 
9j coefficient from the two-nucleon overlap with the ap- 
propriate jj-coupled TNA to construct a set of LS- 
coupled TNA, <£, as 

{ tx S J! ) 

eZ s =hhLS { £ 2 s 32 ) ci T , (4) 

[LSI) 

that satisfy the sum rule 

T,K t ls] 2 = K t ] 2 . (5) 

LS 

Antisymmetry requires, for configurations a where the 
nucleons originate from the same orbital, the [n£j] 2 cases, 



where "F; and <I>( F ) were defined above. The signed 
two-nucleon amplitudes (TNA) C^ T =ci' J//T,T/T will 
be taken from shell-model calculations. They express 
the parentage (amplitudes) for finding each two-nucleon 
configuration a and residue final state / in the over- 
lap with the projectile initial state i, assumed to be 
the ground state. The two-nucleon configurations index, 
a = [fti, 02\ , denotes the spherical quantum numbers of 
the single-particle states occupied by the nucleon pair. 
Hence, /? = ntj. Note that the amplitudes C^ T refer to a 
specific i f, initial to final state transition, which will 
be understood implicitly. 



The details of the shell-model calculations (e.g. the 
model spaces and interactions) used to construct these 
overlaps will be presented with each relevant example in 
the later Sections. 



Expressed in LS'-coupling, the antisymmetrized two- 
nucleon wave function in Eq. (TTJ is 



(2) 



that L+S+T is odd. For nucleons originating from differ- 
ent orbitals this is no longer the case; though for two nu- 
cleons from spin-orbit partner orbitals, [n£j < , nfj>] with 
very similar radial wave functions, the L + S + T =even 
amplitudes are also expected to be significantly sup- 
pressed. 



B. Eikonal model of two-nucleon stripping 

As was developed previously [H, 111 [TTJ] , we exploit 
eikonal reaction dynamics. The elastic ^-matrices de- 
scribing the absorptive interactions of the A-body core 
(in state /) and the two nucleons with the target arc 
calculated in the optical limit of Glauber's multiple scat- 
tering theory [l6|, [13] assuming that these projectile con- 
stituents travel on straight line paths in the interaction 
field of the light target. The reaction is assumed sudden, 
such that the projectile internal co-ordinates are frozen 
on the timescale of this passing and interaction with the 
target. The eikonal S'-matrices are calculated from the 



[^(1)0^(2)]^ =Djij 2 J2 (4Ai4A 2 |LA)(LAS£|/m)L5x ss (1 ! 2)xtt(1,2) 
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nucleoli- and heavy residue-target interactions. These 
interactions were obtained by double-folding the residue 
(core), nucleon (J- function) and the target point particle 
densities with the usual effective nucleon-nucleon inter- 
action, as used elsewhere: e.g. §]■ 

Following Refs. [1, 0, [ll|, from the total absorption 
cross section for the projectile-target system, 



Cabs 



J i Mi J 



(6) 



that includes all events where one or more of the pro- 
jectile constituents are absorbed by the target, we can 
identify and extract the two-nucleon stripping cross sec- 
tion terms, 

ct ^ = -^E /*(*iii5/i 2 (i-i5ii 2 )(i-i5 2 | 2 )i*i) . 

J i Mi J 

n (7) 

As has been discussed elsewhere [llj, the two-nucleon 
stripping probability O str (b, b u b 2 ) = |5/| 2 (1 - |5i| 2 )(l - 
I1S2I 2 ) weights the impact parameters that contribute to 
these stripping events. Stripping requires an absorptive 
(inelastic) interaction of two nucleons with the target, 
but a non-absorptive (elastic) or non-interaction of the 
heavy residue with the target, and strongly localizes the 
reaction to grazing collisions at the projectile surface. 
This simplifies our picture of the reaction mechanism 
and of that part of the overlap function that is probed 
in the knockout reaction. We also contrast the present 
strong surface localisation of the reaction with the term 
peripheral used by some authors to mean impact param- 
eters that sample only the extreme tail (the Whittakcr 
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FIG. 1: Angular momentum and isospin coupling schemes 
used in the calculations. The projectile initial (ground) state 
and final residue state have spin Ji (projection Mi) and J/ 
(projection Mf) and isospins Ti (projection n) and Tj (pro- 
jection Tf), respectively. Each two nucleon configuration a 
involves of a pair of active orbitals j3i and 02 with spherical 
shell-model quantum numbers m, it and ji (projections A; 
and m,i). In LS'-coupling, the nucleon orbital angular mo- 
menta £1 and I2 are coupled to L (projection A), the intrinsic 
spins to S (projection E) and L and S to a total angular mo- 
mentum I, which must also couple the initial and final state 
total angular momenta. 



or Hankel function asymptotics) of the nucleon bound 
state wave function. This is certainly not the case for 
those impact parameters selected by O s tr(b, 61, 62) at the 
projectile energies and with the corresponding absorptive 
S-matrices of interest here. 

Two further assumptions are made. The most im- 
portant, but which is supported by the stripping mech- 
anism's selection of non- absorptive or non-interactive 
events of the residue and target, is to assume there is 
no dynamical excitation/change of state of the reaction 
residue by Sf in the collision - previously termed the 
spectator core approximation. This being the case, 



{& F 'HA)\\S f \ 2 \& F HA)) 



\s r \ 2 s 



FF' 



(8) 



where the bra and ket integrate out the internal coor- 
dinates of the residue and S c is taken to be the residue 
ground state-target elastic scattering S'-matrix. A lesser 
assumption is the heavy-core (no-recoil) approximation, 
that in the required integrals, such as Eq. ([7]), the residue 
impact parameter b c , entering Sf can be replaced by that 
of the center of mass of the projectile b, i.e. b c rj b. 
For the stripping terms this approximation is not in fact 
needed, as a change of integration variable makes it un- 
necessary. 

The result of these assumptions, together with the 
parentage expansion for the two nucleon structure over- 
lap, Eq. (TTJ, is that the exclusive stripping cross section 
to a given final state / can be written 

*£= [db\s c \>± y: <*ri(i-i* 1 2 ) 

J J i MiM f 

(l-|S a | 2 )|*f>) . (9) 

The assumption that the nucleon S'-matrices are spin- 
independent allows one to carry out all spin coordinate 
sums, in preparation for which we separate explicitly the 
nucleon position and spin variable integrations, as 



(F) 



= JdnJ dr 2 (*< F >|...|*f>) 



(10) 



where the final bra-ket term denotes the spin integration. 

Consideration of this (momentum-integrated) strip- 
ping cross section using LS-coupling was made in Ref. 
|8[ , there with an emphasis on the reaction mechanism's 
lack of selectivity in the total spin S of the two nucle- 
ons. In our previous analysis of the longitudinal momen- 
tum distributions (llj the only angular momentum pro- 
jections not able to be summed over algebraically were 
those of the orbital angular momenta of the two nucleons. 

This observation made, we now consider the residue 
longitudinal momentum distributions in the LS- 
representation. The derivation follows a similar pattern 
to that in the jj-coupled algebra and begins from the LS- 
couplcd, spin-integrated, modulus squared of the two- 
nucleon overlap, averaged over initial projections Mj and 
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summed over final projections Mf. Explicitly. 
1 



°i Mi M f 



1 



E c, 

■xcc'TT' 



ITfjl'T' 



(IfiJfMflJiMi) (I'fx'JfMflJiMi) 



MiM f 

(TTT f Tf\TiTi) (T'r'TfTflTiTi 



x([V^(l)®^<(2) 



,T'r' 

I/'m' 



[^(1)8^(2)]?;) . 

(11) 



On performing the sums over Mi and the expres- 
sion is clearly incoherent in the coupled two-nucleon total 
angular momentum J/i, a consequence of the spectator- 
core approximation, Eq. ((HJ). Using the antisymmet- 
ric two-nucleon LS'-coupled forms of Eqs. ^ and (01, 
and assuming the nucleon S'-matrices are also isosphi- 
independent, we can perform the isospin sums with the 
result that Eq. (|TT|) is also incoherent with respect to 
both S and T. Finally, summing over the projections of 
/ and S we obtain the result, incoherent also in L and A, 
namely 



i E 

°i MiM f 



■■^2(TrT f T f \Tinf J2 

T ILSaa 

E, 

A2 A^, 



L 2 



(l 1 \ 1 £ 2 \2\LA)(t 1 \[£' 2 \ l 2 \LA) 



K^O^M - (-) a+ ' J '^i(r 2 )^(ri)] • ( 



fl2l 



The exclusive two-nucleon stripping cross section is then 
given by use of this structure overlap information in Eq. 

In the following we derive explicit expressions for 
the associated exclusive momentum distributions in this 
LS'-representation. 



C. Residue momentum distributions 

Structurally, the expressions for the residue momen- 
tum distributions in the L^-coupling scheme are similar 
to those using jj-coupling The derivations also fol- 
low a largely parallel procedure. The coordinate system 
used is reproduced in Fig. [5] for clarity of the following 
expressions. 

The reaction samples the momentum content of the 
bound-state wave functions of the stripped nucleons in 
the direction of the projectile beam k (i.e. the z-axis). 
For fixed values of the s,, and hence fixed nucleon impact 
parameters bi = \b + Sj|, i = 1,2, this information is 



carried by the functions 



nii) 



dziUp(ri) P\ '(cos0j) exp[iKiZi] , 

(13) 



where is the z-component of the momentum of nucleon 
i in the projectile's rest-frame. We note that the notation 
for TZp(i) is changed from that of Ref. [ll| consistent with 
the notation used for the up{ri). The correct weighting 
of the nucleon absorption probability with the azimuthal 
angle (fi, of Si, is carried by the functions T-L\\'{i) = 
W\\'(b,Si), given by 

Ux\i{i)= / dipi(l - \Si(\b + Si\)\ 2 ) exp[i(fi(X — A')] , 
Jo 

(14) 

with Si the components of the r, in the impact parameter 
plane, i.e. j*j = s.i + z$fc. The remaining details of the 
derivation are completely analogous to those in Ref. [ll[ , 
to which the reader is referred. 

We obtain the projectile rest frame, stripping mecha- 
nism momentum distribution as the incoherent LS and 
isospin decomposition 



da 



if) 



E 



da 



(/) 

LST 



LST 



OrflT (fIT n n 



LSIaa' 



L 2 



J dni J dn2 5(k c + ki+K2) J db |6> c (6)| 2 
E (^iAi^A 2 |LA)(£' 1 A' 1 4A^|LA) 

AAi A2 A^ A2 

J ds\S\ J ds2S2 [direct ~ exchange] , (15) 




FIG. 2: Schematic of the particle coordinates used. Vec- 
tors Si are the components, in the plane perpendicular to 
the beam direction (the z-axis), of the position vectors n of 
the knocked-out nucleons relative to the core of nucleons to 
which they are initially bound. The two nucleons have impact 
parameters bi = b + Si relative to the target nucleus. 
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where the direct term is 

direct = {n XlK (1)7^(1)7^(1)* 

xM W ;(2)^(2)K;?(2)'}, (16) 

and the exchange term, 

exchange =(-l) S+T {Ha 2 a; (1)^(1) 

xH w; (2)^(2)k;?(2)'} . (17) 

It should be noted that Eq. (TT5j) (that contains a fac- 
tor of 2) and these simplified forms for the direct and 
exchange terms compared to [TT| . assume that both the 
integrals over the Ki and the Si will be carried out, and 
so one is computing quantities that are completely sym- 
metric in the two nucleon coordinates. Also, unlike for 
the jj-couplcd scheme, no further re-coupling is required 
to reduce the angular momentum algebra. 

Physically, Eq. (fT5|) shows that the sums of the direct 
and exchange terms over the A arc independent of /. The 
resulting momentum distributions thus depend explicitly 
on L (and S + T via the phase of the exchange term), 
but not on /. We see that the significance of / and the 
nucleon total angular momenta ji is that they will deter- 
mine the relative strengths of the different L and S via 
the amplitudes Thus it is L, and to a lesser extent 

S and T, that will determine the shape of the residue's 
momentum distribution. / on the other hand will be im- 
portant in determining the relative strengths of the L 
and S that contribute. 



III. TWO-NUCLEON CORRELATIONS 

We observe that the expression for residue momen- 
tum distributions is somewhat simpler when using LS- 
coupling, having a more transparent angular momentum 
dependence. However, the dependence on the two nu- 
cleon configurations, via a (and a'), is still less than 
transparent in Eq. (TT5|) . We attempt to elucidate this 
important nuclear structure sensitivity by carrying out 
the A projections sums. Before doing so we introduce 
and discuss the two-nucleon joint position probability 
that summarizes both the strength and the spatial lo- 
calisation (and correlation) of the two nucleons in the 
structure overlaps that affects the stripping yield. 

A. Two-nucleon joint position probability 

We consider the two-nucleon joint position probabil- 
ity relevant to the removal reaction/transition to a given 
final state /, i.e. 

Pf(r 1 ,r 2 ) = ± J2 (^WX ■ (18) 

J i MiM f 



While the production of a given residue final state by 
the two-nucleon knockout mechanism will depend on the 
details of Pf{r\, r 2 ), specifically the extent to which there 
is a spatial proximity of the two nucleons at the projectile 
surface, its overall normalisation and the LS'-composition 
of this normalisation 

Nf = J dr x J dr 2 Pf(r 1 ,r 2 ) 

alT 

LS (alT ) LS 

(19) 

are measures of the likely transition strength. In the 
case of a single (dominant) two-nucleon structure config- 
uration this LS breakdown can also guide the relative 
strengths expected from the different contributing LS 
terms for a given final state. However, when configura- 
tions are mixed or where the initial and final states have 
different parity, interference effects may strongly affect 
these relative strengths. 

Since the projectile is assumed to traverse a straight 
line path in the z-direction, it is useful for what follows, 
and also highly intuitive, to construct the projection of 
the two-nucleon joint position probability onto the im- 
pact parameter plane - the plane perpendicular to the 
beam direction - by integration over the z, of the two 
nucleons, 

Vf(sx,8 2 ) = j "dz 1 J dz 2 p f (r 1 ,r 2 ) . (20) 

The relevant spatial correlation for the reaction is now 
the degree of localisation of the probability with respect 
to the two nucleon coordinate projections s, in this im- 
pact parameter plane. 

In what follows the correlation of the two nucleons is 
concisely expressed as a function of the angular separa- 
tion, u, of their position coordinates r^. Clearly, the Zi 
integrated joint probability Vf(s\, So) will see a smeared 
version of this correlation function since fixed Si will sam- 
ple a range of ui. However, since the reaction is sur- 
face localized and the target is light (small) the effective 
thickness in the Zi will tend to be rather restricted and 
Vf(si,s 2 ) will remain a useful construct and intuitive 
link to the magnitudes of the two-nucleon knockout cross 
sections. 

As was indicated by Figs. 2 and 3 of Ref. [l8[, and 
will be emphasized here, the total angular momentum of 
the final state and the detailed TNA of the wave function 
can strongly affect the two-nucleon joint position prob- 
ability, its projection and the magnitude of the removal 
cross sections. The (shell- model) structural correlations 
may also enhance or suppress particular total orbital an- 
gular momenta and so may affect the residue momentum 
distributions also. 
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B. Angular Correlations 

Despite the relative simplifications introduced by LS- 
coupling, the momentum distribution expression, Eq. 
(|15p . remains a complicated weighted sum of wave func- 
tion transforms. Moreover, it still depends on the orbital 
angular momentum projections. To clarify the underly- 
ing sensitivity to two-nucleon correlations we simplify the 
spin-integrated two-nucleon joint position probability of 
Eq. (fl"8)) by summing out the A projection labels. 

The relevant terms we need to simplify are, for the 
direct terms of Eq. (fT2j) , 

r^ 2 ^(fi,r 2 )= £ (W 2 A 3 |LA)(4Ai4Ai|£A) 

AAi A 2 A'j 

(21) 

Combining the spherical harmonics of the same argu- 
ment, summing the A projections, and using the spherical 
harmonics addition theorem one obtains 

x(-l) fe (£ 1 0£' 1 0|fcO)(€ 2 0£20|fcO)P fc (cosw), 

(22) 

where w is the angular separation of the two nucleons, 
i.e. cos uj = T\ ■ r 2 /r\r 2 . A similar result can be found in 
Ref. [Uj]. 

The angular correlation function ^f g g/ a (oj) is seen 
to be independent of the total angular momentum / and 
of the individual angular momenta ji of the nucleons. 
However, it depends explicitly on their orbital angular 
momenta and on the total orbital angular momentum L. 
The form written above is that for the direct terms of 
Eq. (|T2|) . The exchange terms differ by a phase due to 
the reordering of the angular momentum labels in the 
exchange form of Eq. (f2~Tj) , as is given below. 

The radial behaviors associated with the direct and 
exchange terms of the joint-probability density are 

Ua a ,(ri,r 2 ) = u^(ri) u 02 {r 2 ) u P[ (n) u^(r 2 ) 
+ up 2 (r 1 ) u Pl (r 2 ) u 02 (n) u^(r 2 ) , 

Ua a >( r i,r 2 ) =up 1 (ri) u fj2 (r 2 ) w^(n) u#(r 2 ) 

+ W/3 2 (ri ) Ufa (r 2 ) (n ) up> 2 (r 2 ) . (23) 

In terms of these and the corresponding direct and ex- 
change angular correlation functions, the two-nucleon 
joint-probability density is 

^ ^ <T /T <T /T D D , 
Pf(rx,r 2 ) = Y, £ 5 f 2 a a {TrT f r s \T in f 

LST Iaa' ^ 

x[u° a ,( ri ,r 2 )r L - D (uj) 

-{-) S+T U* a ,{n,r 2 )T L > E {u)], (24) 



with T l - d (uj) = T% itae , given by Eq. fl22J) and 

T L ' E {u) = {-Y'i +t '*- L T^ w {u). (25) 

It is clear therefore that the angular correlation func- 
tion dictates how the spatial correlations change with 
angular momentum coupling, and that L is crucial, the 
U{ri,r 2 ) being dependent on the Pi but independent of 
the angular momentum coupling. Clear also is that, in 
cases where the radial wave functions for all active orbits 
are similar, the angular correlation function alone will 
determine the differences in residue momentum distribu- 
tions for the different possible angular momentum cou- 
plings. As was discussed earlier, these differences, gen- 
erated at the angular correlation function and the two- 
nucleon density level, will be more distinct than in the 
projected density, Eq. (|2T)|) where fixed co-ordinate pairs 
(si,s 2 ) sample a range of angular separations uj and so 
will smear the spatial correlations. 

Uncorrelated two nucleon models, discussed in the In- 
troduction and Refs. [1, [8|, that neglect antisymmetriza- 
tion, angular momentum coupling and parentage coeffi- 
cients lead to a constant, w-independent correlation func- 
tion. For two nucleon removal from a single [si/2] 2 con- 
figuration, the angular correlation function is also seen 
to be w-independent (k = 0) and the uncorrelated (see 
Ref. (Tlj ) and fully correlated residue momentum distri- 
butions will be identical. 



C. Cross shell excitations 

Here we consider briefly the implications for two- 
nucleon knockout from configurations with j3\ and fi 2 of 
different parity. It is well established that the addition of 
shell- model configurations with lhio, ituo single particle 
excitations are required to obtain a high degree of surface 
pairing (see e.g. |20l - [25| ). 

We obtain a similar result here, by considering the 
symmetry of the angular correlation function about u> = 
ir/2. In r^^, e , (uj), only the Legendre polynomial de- 
pends on uj, with the property that Pfc(cos[7r — uj]) = 
(— l) fe Pfc(cos uj). Since the values of k are restricted to 
be odd or even by the parity Clebsch Gordan coeffi- 
cients, the angular correlation will be even about n/2 
for 7T£7T£' = +1 and odd about 7r/2 for t:(t:i' — —1. In 
the absence of single-particle excitations of the kind \ fvui, 
the probability for finding the nucleon pair with angular 
separation ui = and 7r are equal and a high degree of 
two-nucleon pair/cluster structure will not be obtained. 

So, pair correlations will be enhanced in cases when 
there is mixing between two-nucleon configurations 
where the orbital angular momenta are of different par- 
ity. Whether the interference is constructive or destruc- 
tive will depend on the sign of ^f lg2£ ' e i (w) near u) = 0, 
the relative signs of the 9j coefficients, and the relative 
signs of the TNA. A specific two configuration example 
will be presented in the Section IIV CI below. 
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These results are quite general in that they do not 
depend on the pair total angular momentum /; enhance- 
ments in the spatial correlations in the two-nucleon den- 
sity may be found for I 7^ 0. 



IV. ILLUSTRATIVE EXAMPLES 



TABLE I: Isospin format two nucleon amplitudes for the first 
and second lc, B(f + ) T = states populated by neutron- 
proton removal from 12 C. The relative strengths of the two 
nucleon amplitudes are different for the two states. 



T 7 ' 



[1P3/2 



[lp 



1/2] 



0.69899 
-1.13385 



0.97868 
0.22886 



-0.01067 
0.36314 



Previous calculations of exclusive two-nucleon removal 
residue momentum distributions noted a strong sensitiv- 
ity to the total angular momentum of the removed nu- 
cleon pair. Here, by writing this momentum-differential 
cross section in LS- coupling, and by a consideration of 
the angular correlations inherent in the two-nucleon joint 
probability function, it becomes apparent that the cru- 
cial sensitivity of this observable is to the total orbital 
angular momentum values, L, contributing to the transi- 
tion. These different L components will contribute inco- 
herently to the cross section yields and their momentum 
distributions. These theoretical observations and the re- 
sulting sensitivity of the momentum distribution observ- 
able offers the potential to probe more subtle features 
of the nucleon pair's configurations and the correlations 
present in the shell-model wave functions used. 

A generic first example will arise if the predominant 
two-nucleon configuration populating a given final state 
involves one of the nucleons in an s-wave orbital. In this 
case the total orbital angular momentum is restricted to 
the orbital angular momentum of the second active or- 
bit, L = £, and thus L is pure. It is expected therefore 
that there can be distinct differences in the momentum 
distributions, even for states of the same total angular 
momentum I. For example, for two 3 + final states built 
from [<77/2, S1/2] an( l from [^5/2,51/2]- More generally, 
even where there is significant mixing and several active 
configurations, the structure of specific states in the spec- 
trum can be rather L-pure. So, the reaction will proceed 
by a particular L with a momentum distribution that is 
characteristic of this structure. 

In the following we discuss specific examples from dif- 
ferent A and Z regions of the nuclear chart. In each 
example the nucleon bound state radial wave functions 
required for the two-nucleon overlaps are calculated us- 
ing a Woods-Saxon potential well with a spin-orbit term 
of depth 6 MeV and a diffuseness parameter ao=0.7 fm. 
Unless stated otherwise, the geometries (the radius pa- 
rameters ro) of the potential wells in each case were ad- 
justed to reproduce the root mean square radii and the 
separation energies of spherical Hartree-Fock calculations 
using the Skyrme (SkX) interaction parameterization (26j 
for the active orbitals in question. The specific procedure 
was detailed in Ref. Q . These fitted geometries are then 
used to calculate the radial wave functions needed us- 
ing the empirical, effective nucleon separation energies. 
Where required, shell-model calculations are performed 
using the code OXBASH [27} • The model spaces and in- 
teractions used are specified for each case studied, below. 



TABLE II: Two-neutron stripping cross sections gls for pop- 
ulating the first two T = 0, J/ = 1 + final states in 10 B for a 
12 C beam energy of 2100 MeV per nucleon. All cross sections 
are in mb. 



Jf 


(701 


<T10 


Oil 


(721 




It 


2.41 


0.00 


0.00 


0.06 


2.47 


It 


0.60 


0.59 


0.00 


0.63 


1.81 



A. p-shell example: 12 C(-np) 

Here we consider the removal of a (T = 0, 1) neutron 
and proton (np) pair from 12 C at 2100 MeV/nucleon on 
a 12 C target. The proton and neutron orbits are taken 
to be identical with radial wave functions calculated in a 
Woods-Saxon potential, using an average nucleon charge 
Z = 0.5. The geometry of the Woods-Saxon potential 
was fixed with ro = 1.31 fm, ao = 0.7 fm. Both the 
10 B residue and 12 C target were assumed have Gaussian 
shaped mass distributions, with rms radii 2.30 and 2.32 
fm respectively. The isospin format TNA are calculated 
using OXBASH in a p-shell model space using the wbp 
interaction [28| , as in previous studies 0, |3(| • A more 
complete consideration of two-nucleon removal from 12 C 



will be discussed in a forthcoming paper [31 1 . 

As a specific example, we consider the first and second 
T = 0, 10 B(1+) final states. The TNA for these states 
are shown in Table |TJ The relative magnitudes of the 
contributing two-nucleon configurations to these states 
are different and it is of interest to consider how these 
differences might affect the cross sections and their mo- 
mentum distributions. The sum of the squared TNA for 
the first and second states are 1.45 and 1.47, respectively, 
thus in the absence of interference terms the incoherent 
sum of contributions from each of these configurations 
would yield very similar cross sections. 

That this is not the case is shown by the calculated 
two-nucleon stripping cross sections presented in TablclHl 
The calculated momentum distributions are also rather 
different, as is shown in Fig. [3] 

These differences can be understood by reference to 
the projected two-particle joint position probabilities for 
the two states, which are strikingly different. The first 
1 + state shows strong spatial localisation of the two- 
nucleons, favorable for the two-nucleon removal cross sec- 
tion. Both example position probabilities manifest the 
expected symmetry about a nucleon angular separation 
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0.6- 



present model space, independent of the configurations 
(a, a') of the pair. It follows that the weight of each 7 
and S term in a state of given 7 and T is proportional to 



P 



IT 
LS 



IT rflT 



LS^a'LS 



D a D a , 



(26) 



0.4 



T3 



0.2- 
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FIG. 3: Normalized residue momentum distributions for the 
first (solid) and second (dashed) 10 B( Jf = l + ) states populated 
m np knockout from 12 C at 2100 MeV per nucleon. Though 
the same two-nucleon configurations contribute to each state, 
the differently weighted TNA result in distinct momentum 
distributions. 




-4 -2 2 4 
s x (fm) 



-2 2 
s x (fm) 



FIG. 4: (Color online) Impact parameter plane-projected 
joint position probabilities for (a) the first and (b) second 
T — 10 B(1 + ) states populated via np knockout from 12 C. 
The plot shows the impact parameter plane probability dis- 
tribution of nucleon 2 for nucleon 1 positioned at s x = 2.9 
fm, s y — fm. The spatial correlations of the nucleon pairs 
in these two states are fundamentally different, leading to 
markedly different momentum distributions; see Fig. [3] The 
color scale (white-blue-green-yellow-red-black) is common to 
both plots. 



of = 7r/2, since the model space is restricted to the 
p-shell and the active orbitals have the same parity. 

We can extend this p-shell example further to illustrate 
the potential for large sensitivity to the underlying struc- 
ture. It is clear from Eq. (f2~4")l that within a p-shell model 
space the relative strengths of different LS combinations 
are determined solely by the TNA and the nucleon config- 
urations involved. So, neglecting any minor differences in 
the p-wave radial wave functions, due to spin-orbit split- 
ting, the entire square bracketed term in Eq. (|24[) is in- 
dependent of the total angular momenta {ji} and, in the 



In the sprit of studying the extremes of possible sensi- 
tivity of the momentum distributions, we may force any 
one of these P/g to be zero, and solve for the relative 
strengths and phases of the £^,"£5 and C^ T needed to 
achieve this. 

Fig. [5] illustrates such examples for assumed 7=1, 
T = states populated via the configurations o,\ — 
[lpi/2] 2 and <X2 = [l.P3/2j lfi/2]- Calculations for two sets 
of TNA are shown; one set chosen to eliminate L = 0, S = 
1 contributions (requiring = —4(7^°, dashed curve) 
and the other to eliminate L = 2, S = 1 contributions 
(requiring = C^/2, solid curve). These different 1 + 
state TNA produce wide and narrow residue momentum 
distributions, respectively, the difference in the FWHM 
widths being almost a factor of two. The figure also 
shows the 7 = 2, T = 1 momentum distribution (open cir- 
cles), populated via a 2 = [IP3/2, lpi/2] and a 3 = [lp 3 / 2 ] 2 , 
where the TNA were chosen to eliminate L = 2,5= 1 
contributions (requiring = — \/2Cal)- Once again 
this gives a relatively narrow distribution and, moreover, 
this I ~ 2 distribution is narrower than that for the 
(7 = excluded) 7=1 distribution (dashed curve) de- 
scribed above. These examples break the tie between the 
width of the momentum distribution and the 7 value of 
the transferred pair. Whether or not nuclear states with 
these TNA are physically realized, these limiting cases 
demonstrate how details of the microscopic structure of a 
given state may strongly influence the shapes and widths 
of the expected residue momentum distributions. 

Similarly, we note the expectation that transitions to 
7 = 2, T = and 7 = 3,T = states of 10 B will, within 
a p-shell model space, yield identical theoretical momen- 
tum distributions since both transitions are pure 7 = 2 in 
nature. In this instance also the width of the momentum 
distribution does not provide a direct measure of Jf = 7. 

We note that consideration has only been given to the 
direct population of 10 B. In principle, indirect popula- 
tion by single nucleon knockout followed by evaporation 
of the unlike nucleon may be possible, although we ex- 
pect this indirect, two-step pathway to be very weak, due 
to the large nucleon separation energies in the relevant 
A = 11 systems and the very small predicted shell model 
strength for one nucleon removal to states above these 
first A = 11 nucleon thresholds. High precision (stable 
beam) observations of final state exclusive 10 B momen- 
tum distributions would clarify such aspects of the reac- 
tion mechanism that are currently assumed. 
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K c (MeV/c) 

FIG. 5: Theoretical 10 B residue momentum distributions at 
2100 MeV per nucleon. Shown are the expectations for two 
/ = 1 T = states where the TNA have been tailored to 
exclude L = 2, S — 1 (solid line) and L — 0, S — 1 (dashed 
line) contributions from the two-nucleon wave function, re- 
spectively. The line with open circles shows the expectation 
for an I = 2 T = 1 final state where the TNA were similarly 
chosen to eliminate L = 2, S — 1 contributions. 



B. sd-shell: 28 Mg(-2p) and 26 Si(-2n) 

Exotic nuclei with valence nucleons in the sd-shell have 
been the focus of several two-nucleon removal experi- 
ments, studying the evolution of structure away from the 
valley of /3-stability. The initial and final state struc- 
tures are often well described within conventional sd-shell 
model space calculations offering good test cases for stud- 
ies of the reaction mechanism. Details of their residue 
momentum distributions could offer an additional test 
of the shell model and the reaction mechanism in this 
region. 

We first review the two-proton knockout from 28 Mg at 
83.2 MeV/nucleon on a 9 Be target, previously studied 
in Refs. H, H, [ll|. To date, this is the only experi- 
mental example with measured final-state exclusive 26 Ne 
momentum distributions. Four states were populated, 
being the + ground state, the first and second 2+ states 
and the first 4 + state. Previous work demonstrated the 
significant difference between the ground state and 4 + 
residue momentum distributions, despite strong experi- 
mental (reaction target) broadening of the measured dis- 
tributions. 

We comment here on the effects on the 2.02 and 3.70 
MeV 2 + state momentum distributions of the subtle dif- 
ferences in their TNA, tabulated in Rcf. ||. To remove 
the small difference in the average separation energies of 
the protons for the two states, calculations used iden- 
tical radial wave functions, but this binding effect is in 
practice negligible. The calculated widths of the residue 
momentum distributions arc different by ~ 10%. Clearly 



a higher statistics experiment would be required to ex- 
amine this difference predicted by the shell model. There 
are however other examples where the sii-shell model pre- 
dicts TNA that exhibit a larger degree of sensitivity, as 
e.g. the following. 

A second specific example is the two- neutron (T = 1) 
knockout from 26 Si, measurements for which were re- 
ported in Ref. [HH; made at 109 MeV per nucleon on 
a 9 Be target. Details of the nucleon radial wave func- 
tions and S'-matrices can be found in Ref. @. Popu- 
lations of two excited states in 24 Si were observed, the 
first 2+ state at 1.86 MeV, and a state at 3.41 MeV cor- 
responding to a theoretically-predicted (2 + , 4+) doublet, 
with theoretical excitation energies of 3.867 and 3.962 
MeV. The cross sections for these measured and theo- 
retical states were analyzed Q assuming that the second 
excited state was the second 2 + state. Momentum dis- 
tributions, if available, would easily distinguish between 
such 1 = 2 and 1 = 4 possibilities. Our interest here 
is more subtle. We consider the expected differences in 
the momentum distributions of the two 2 + states arising 
from their underlying sd-shell model structures. 

The TNA were calculated using OXBASH within an sd- 
shell model space using the USD interaction [33| and 
are presented in Table IIIII The TNA calculated using 
the USDA and USDB interactions j34| were found to 
be very similar to the USD values. Both states have 
mixed sd-shell configurations. Inspection of the TNA 
might suggest that since the second state has a stronger 
[2s!/2, ld.5/2] configuration it may favor L = 2 more 
strongly, but there is significant mixing. 

Despite the strong mixing in both states, the shell- 
model TNA predict each state to be populated predom- 
inately by a single and distinct total orbital angular mo- 
mentum L, L = 1 and L = 2, respectively. The cal- 
culated LS'-coupled two-nucleon stripping partial cross 
sections reveal this, as are shown in Table HVl The dom- 
inance of L = 1 and L = 2 in these states generates 
the significantly different 2 + state momentum distribu- 
tions of Fig. |51 the l\ state having a 30% larger width. 
Exclusive measurements for these states would not only 
clarify if the second excited state is the 2^, but could 
also confirm the L = 2 dominance prediction of the sd 
shell-model calculations. 

To consolidate our understanding of such sensitivity, 
we consider a further simplified example where a sin- 
gle configuration is expected to dominate. We consider 
the two configurations [ld 5 / 2 ] 2 and [lds/2, 1^3/2], both of 
which can contribute to 4 + states. States with such sim- 
ple configurations may not be realized in 2 Si, since 4 + 
states in 24 Si are thought to be unbound, but the exam- 
ple will serve to illustrate the expected differences that 
may occur elsewhere in the sd-shell. 

We construct the TNA as C 41 = such 
that Nf = 1, see Eq. (fl"9)) , and the resulting LS- 
decomposition of strengths is given in Table [V] It is 
very clear that the [lds/2, ld 3 / 2 ] configuration weights 
L = 4 significantly more strongly than does [ld 5 / 2 ] 2 
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TABLE III: USD shell-model TNA for the first two 2+ states in 24 Si, populated in two-neutron knockout from 26 Si. 



Jf 


[W5/2] 2 


[ld 5 / 2 , ld 3/2 ] 


[W3/2] 2 


[2s 1/2 , ld 3/2 ] 


[2si /2 , ld 6/2 ] 




-0.70074 


0.43499 


0.00594 


-0.00188 


-0.02781 


2i 


-0.38021 


-0.12354 


-0.12945 


-0.15876 


-0.58292 



TABLE IV: Two- neutron stripping cross sections, ijls, f° r 
population of the first two 2 + final states in 24 Si. All cross 
sections are in mb. 



Jf 




O"20 


C"21 


0-3 1 


<?str 




0.17 


0.02 


0.00 


0.00 


0.19 




0.01 


0.17 


0.01 


0.00 


0.19 




k c (MeV/c) 

FIG. 6: Normalized residue momentum distributions for the 
first (solid) and second (dashed) 24 Si(2 + ) states populated 
in 2n knockout from 26 Si. Though the same sd-shell two- 
nucleon configurations contribute to each state, their TNA re- 
sult in distinct L makeup and momentum distributions. The 
(FWHM) peak widths are are 250 and 330 MeV/c, respec- 
tively. 

and the expectation is a wider momentum distribution. 
As noted in Section III Al in this case we would expect 
L + S + T=evcn contributions to be significantly sup- 
pressed due to the two-neutron antisymmetry and the 
similarity of the radial wave functions for the active spin- 
orbit partner orbitals. This is indeed the demon- 
strated by the stripping cross sections of Table [V] The 
estimated strengths, Nj S , are seen to give a reasonable 
guide to the expected cross sections for these single con- 
figuration examples. 

The results of the calculations, shown in Fig. [7J con- 
firm the differences in the momentum distributions ex- 
pected from our simple consideration of the Nj S . Again, 
the specifics of the underlying structure predict consid- 
erable and observable differences in the expected residue 
momentum distributions. 



TABLE V: LS-coupled strengths, Nf s , and LS partial strip- 
ping cross sections for pure [lds/2] 2 and pure [lds/2, ld 3 /2] 
two- neutron configurations populating 4 + final states in 24 Si. 
The cross section for L = 4, S — 1 is negligible, see text. 



Q 


ATf 


Nf 


Nf 


031 


040 


(741 


[lds/2] 2 


0.8 


0.2 


0.0 


0.23 


0.09 


0.00 


[Ids/a] [Ids/a] 


0.1 


0.4 


0.5 


0.06 


0.35 


0.00 




k c (MeV/c) 

FIG. 7: Examples of the projectile rest frame residue momen- 
tum distributions for I — A states in the 26 Si(— 2n) reaction, 
arising from different neutron pair configurations. The solid 
lines and open points assume [lds/2] 2 neutron removal, the 
dashed lines and solid points assume [ld 5 / 2 , ld 3 / 2 ] neutron 
removal. The full distributions are normalized to 1 at « c = 0, 
with each contributing L partial distribution scaled by the 
same factor. Circles show L = 3, S = 1 contributions and 
squares show L — 4, S = 0, with the total shown by the line. 
The different relative strengths of L = 3 and L — 4 for the 
two configurations generates significantly different 24 Si mo- 
mentum distributions for states of the same I (the dashed 
and solid curves). 

C. Cross shell: 54 Ti(-2p) 

This (T = 1) reaction, reported in Ref. [35j . demon- 
strated the potential for two-nucleon knockout to probe 
cross-shell proton excitations in neutron rich nuclei. In 
particular, a 52 Ca(3~,3.9 MeV) state was populated in 
two proton removal from 54 Ti(0 + ) on a 9 Bc target at 72 
MeV per nuclcon. Details of the eikonal S-matrices and 
nucleon radial wave functions can be found in Ref. 13511 . 
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Previous theoretical estimates for the 3 state yield as- 
sumed pure [I/7/2, 2si/ 2 ] or [1/7/2,1^3/2] configurations, 
providing an estimated upper limit for the cross section 
to this state as an incoherent sum of these contributions. 

Taking instead a coherent sum will give (a) a different 
total cross section, and (b) a different residue momentum 
distribution. Here we assess the expected sensitivity to 
the relative strengths and phases of these two configura- 
tions. We calculate the two-proton stripping cross sec- 
tions and momentum distributions as a function of the 
TNA for [l/ 7/2 , ld 3/2 ] removal, Cf d , and for [l/ 7/2 , 2s 1/2 ] 
removal, Cjl . For either of these pure configurations the 
stripping cross sections scale with [CjJ] 2 . To maintain 
an overall scaling when the configurations are mixed, the 
two amplitudes are adjusted such that 

[Cjl? + [Cf d ] 2 = 1 , (27) 

with Cfl assumed positive. The total incoherent strength 
thus remains constant. We calculate the stripping cross 
sections for each contributing LS combination and the 
full-width half-maximum (FWHM) for the residue mo- 
mentum distributions. This is not the whole story for the 
momentum distribution - there are also subtle changes 
of shape beyond the nominal width - but this FWHM 
width provides a guide to the expected behavior. The 
resulting calculations are shown in Fig. |SJ A few points 
follow immediately; the [I/7/2, 2S1/2] configuration only 
contributes to the L = 3 cross section, giving no inter- 
ference with L = 2 and L = 4. So, these latter terms are 
simply proportional to [Cj- d ] 2 and are zero at the centre 
of the plot. The L = 4 contributions are also generally 
weak and the overall width of the residue momentum dis- 
tribution is largely determined by the relative strengths 
of the L = 2 and L = 3 contributions. 

Both the cross section and FWHM of the momentum 
distribution show a strong sensitivity to the mixing of 
the two configurations; the cross section varies by a fac- 
tor of two and the width of the momentum distribution 
by 25%. It is clear that the underlying structure and the 
relative strengths of the two-nucleon amplitudes are crit- 
ical to determining both the removal cross section and 
the shape of the momentum distribution. Of interest are 
the extremes of the plot, with \Cj d \ ss 1. Here both the 
two-nucleon removal cross sections and momentum dis- 
tribution widths are acutely sensitive to the small admix- 
tures of the [I/7/2, 2s!/ 2 ] configuration, but also strongly 
dependent on its sign. If the two amplitudes are of oppo- 
site phase then both the cross section and width decrease 
rapidly. Conversely, they increase rapidly if in phase. 
This is indicative of a sensitivity to small cross-shell ad- 
mixtures in many cases. 

We now contrast the impact parameter plane- 
projected two-nucleon joint probability distributions, for 
values of Cf d = ±y/2/2, in Fig. ED The difference in 
this cross-shell case is now striking and it is clear that 
taking the amplitudes to be in phase (the +ve choice) 
enhances the two-nucleon spatial correlations, that then 




FIG. 8: (a) Two-nucleon stripping cross sections and (b) the 
full-width at half-maximum of the residue momentum distri- 
bution when populating the 52 Ca(3~, 3.9 MeV) state. These 
observables are shown as a function of the amplitude C^\. 
The (positive) amplitude Cj\ of the second configuration is 
given by Eq. 1271 The lower panel shows the partial cross sec- 
tions with {L, S}, for the values {2, 1} (dashed), {3, 0} (solid), 
{3, 1} (dotted) and {4, 1} (dot-dashed). The total stripping 
cross section is shown by the solid line with open triangles). 



drives the significantly larger two-nucleon knockout cross 
section that is calculated and differences in the residue 
momentum distribution. 

As is clear from Fig. [8j a precise measurement of the 
residue momentum distribution for this reaction would 
allow an estimate of the relative strengths and phases of 
the amplitudes of the two assumed active two-nucleon 
(cross-shell) configurations. The knockout cross sections 
themselves are also shown to depend strongly on the mix- 
ing. To date, analyses of two-nucleon knockout from 
exotic (asymmetric) systems have shown that the the- 
oretical cross sections overestimate those measured ex- 
perimentally by of a factor of about two, quantified as 
l^theor ~ 0.55; see e.g. Ref. [1]. This 
suppression effect thus introduces an ambiguity in the ab- 
solute cross sections that is significant at the level of the 
differences being shown in Fig. [8j Such suppressions, of 
the cross sections predicted using the shell-model spec- 
troscopy, may themselves be, at least in part, a man- 
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FIG. 9: (Color online) Impact parameter plane-projected two- 
nucleon joint position probabilities, for L = 3, S — 0, for 
(a) Cf d = s/2/2 and (b) Cf d = -\/2/2 for the 52 Ca(3", 
3.9 MeV) state populated in two-proton knockout form 54 Ti. 
The plot shows the probability distribution of nucleon 2 when 
nucleon 1 is positioned at the back circle. The color scale 
(white-blue-green- yellow-red-black) is the same for both plots. 
The source of the differences in the calculated cross sections 
(see Fig. [HJ for these two choices of TNA is evident in the 
pair correlations manifest in these projected two-particle joint 
position probabilities Note the asymmetry at <pi2 = 90° (i.e. 
s x =0), induced by mixing of different parity orbitals. 

ifestation of the use of TNA calculated in a truncated 
shell-model space and that exclude a large number of 
(small amplitude) cross-shell configurations. Based on 
the limited measurements available to date, there is no 
indication that the (missing) physics that drives the sup- 
pression of cross-section strength has implications for the 
shape of the residue momentum distribution. Additional, 
more accurate exclusive final state data are needed to as- 
sess these expectations further. 

V. SUMMARY 

We have discussed the momentum distributions of the 
heavy residues after two-nucleon knockout reactions us- 
ing L.!?- coupling. The main factor determining the width 
of these momentum distributions is shown to be the 
transferred total orbital angular momenta L of the two 
nucleons. We provide insight into the expected widths of 
momentum distributions from the removal of the nucleon 
pair from different configurations showing that informa- 
tion can be gained from and upon the strengths of the 
theoretical two nucleon amplitudes and the contributing 
L they generate. The unambiguous observation of effects 
associated with specific pairs of nucleon orbitals may re- 
quire transitions to final states that are relatively pure or 
simple configurations. Some illustrative examples were 
presented and discussed. 

The conclusion of previous work - that the residue mo- 
mentum distribution was simply characteristic of the fi- 
nal state spin - has been considered in further detail. It 



is true that, generally, higher spin final states will lead to 
wider residue momentum distributions, but that the de- 
tails of the shell model two-nucleon overlap are important 
in understanding the details of the residue momentum 
distributions. Quantitative testing and confirmation of 
such sensitivity to the underlying structure will be essen- 
tial for the exploitation of two-nucleon knockout methods 
and their extension for deformed nuclei. 

The critical importance of configurations of different 
parity in enhancing pairing correlations is demonstrated 
by consideration of the angular correlations inherent in 
the two-nucleon density. Whilst discussed here in the 
context of two-nucleon removal reactions and enhance- 
ments of two-nucleon removal cross sections, such consid- 
erations, of large basis TNA, and the importance of small 
admixtures of different parity is entirely general. In the 
context of the suppression of shell model strength, pre- 
vious studies with radioactive beams have demonstrated 
that the theoretical cross sections overestimate experi- 
ment by about a factor of two. It will be important to 
experimentally verify the influence of cross-shell excita- 
tions on structurally better-understood cases, such as for 
12 C, 16 and 40 Ca, to clarify the extent to which the nec- 
essary reductions may depend (in part) on the truncated 
model spaces used. It will also be important to further 
assess the importance of cross shell proton-excitations in 
studies of islands of inversion using the two-proton knock- 
out methodology (see e.g. [36l438| ). where very strong 
reductions of two-proton knockout cross sections are ob- 
served. 

Here our emphasis has been on light and medium 
mass projectiles. Another interesting example is the two- 
proton removal reaction from 208 Pb; not only are there 
a large number of active orbitals, producing a plethora 
of states, but the majority of states are good two-proton 
hole configurations with minimal mixing. 

The study of such reactions with odd-mass projec- 
tiles brings an added layer of complication with, typi- 
cally, each final state being populated via several nucleon 
pair total angular momenta. The widths of the residue 
momentum distributions are then no longer simply re- 
lated to a single final state spin. However, the under- 
lying structure sensitivity discussed here may still yield 
characteristic widths for different final states in the same 
residue, somewhat independent of the final state spin. 
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